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Abstract 

A manifestly covariant, coordinate independent reformulation of the Ther- 
modynamic Field Theory (TFT) is presented. The TFT is a covariant field 
theory that describes the evolution of a thermodynamic system, extending the 
near-equilibrium theory established by Prigogine in 1954. We introduce the 
Minimum Dissipation Principle, which is conjectured to apply to any system 
relaxing towards a steady-state. We also derive the thermodynamic field equa- 
tions, which in the case of a— a and (3— 13 processes have already appeared in the 
literature. In more general cases the equations are notably simpler than those 
previously encountered and they are conjectured to hold beyond the weak-field 
regime. Finally we derive the equations that determine the steady-states as 
well as the critical values of the control parameters beyond which a steady- 
state becomes unstable. 
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1 Introduction 



In 1999 the Thermodynamic Field Theory (TFT) was proposed, which describes 
the behavior of thermodynamic systems beyond the hnear regime 0-0. This ap- 
proach starts from the concept of entropy production, which can always be expressed 
as a bihnear form of some variables characterizing the departure from equilibrium. If 
these variables are the thermodynamic forces, it is possible to formulate a covariant 
TFT and to find the thermodynamic field equations whose solutions give the general- 
ized relation between the thermodynamic forces and the conjugate flows. 

In refs 0-0, the thermodynamic field equations have been obtained from three 
postulates: 

1) The shortest path principle 

2) The closedness of the thermodynamic field strength 

3) The principle of least action 

where the thermodynamic field strength is determined from the skew-symmetric part 
of the tensor relating the thermodynamic forces with the conjugate flows. The second 
postulate reflects the observation that, in all cases examined so far, the sources of the 
thermodynamic fields i.e., the internal fluctuations, time-dependent boundary condi- 
tions and external perturbations, only affect the balance equation of the dual form of 
/. This experimental observation allows us to identify the system of thermodynamic 
fleld equations as an analogue of Maxwell's equations. In the present note we will 
see that the second postulate is analogue to the absence of magnetic sources in a 
higher-dimensional version of Maxwell's equations. 

The validity of this theory in the weak-field approximation has been tested in 
many examples of a — a or /? — /3 processes such as, for example, the thermoelectric 
effect and the unimolecular triangular chemical reaction (see refs jl| and 0). In the 
first case, the TFT provides new predictions when the material is subjected to strong 
electric field and, in the second, the thermodynamic field equations reproduce the 
well-known De Bonder law. 

Transport processes in magnetically confined plasmas, which are of a — a type, 
have also been analyzed using the TFT. In particular, in ref. 0] the thermodynamic 
field equations were solved in the weak-field approximation of the classical and the 
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Pfirsch-Schliiter regimes. We found that the TFT does not correct the expressions 
for the ionic heat fluxes predicted by the neoclassical theory. On the other hand, 
the fluxes of matter and electronic energy (heat flow) are enhanced in the nonlinear 
classical and Pfirsch-Schliiter regimes. This phenomenon would have been amplified 
had we used the strong (exact) thermodynamic field equations. These results are in 
line with the experimental observations. 

The thermodjTiamic field equations, in the weak-field approximation, have also 
been applied to several a — (3 processes. For example the Field-Koros-Noyes model, 
in which the thermodynamic forces and flows are related by an asymmetric tensor, 
was analyzed in ref j7|. Even in this case the numerical solutions of the model are in 
agreement with the theoretical predictions of the TFT. Recently, the Hall effect j^, 
jof and, more generally, the galvanomagnetic and thermomagnetic effects have been 
analyzed in nonlinear regimes. In each of these papers it was shown that the TFT 
successfully describes the known physics in the nonlinear regimes and also predicts 
new interesting effects such as the nonlinear Hall effect. The theoretical predictions of 
the nonlinear Hall effect have been confirmed experimentally. In ref ^ it was shown 
that, for materials with low thermoelectric power coefficients and in the temperature 
range of available experimental data, the theoretical predictions of TFT agree with 
experiments. 

When the thermodynamic system approaches equilibrium, the solution of the field 
equations reduces to the Onsager-Casimir tensor and satisfies the Prigogine theorem 
of minimum entropy production 0]. Far from equilibrium Glansdorff and Prigogine 
have demonstrated that, for time independent boundary conditions, a thermodynamic 
system relaxes to a steady-state satisfying the Universal Criterion of Evolution 
Glansdorff and Prigogine obtained this result in 1954 using a purely thermodynamical 
approach. As indicated in refs 0]- ^, making use of the second principle of thermo- 
dynamics, the Shortest Path Principle should ensure the validity of the Universal 
Criterion of Evolution. 

In the present manuscript, using the coordinate independent language of Rieman- 
nian geometry, we reformulate the TFT in a generally covariant way. This allows us 
to extend several weak-field results to the strong-field regime. After expressing the 
concepts of thermodynamic forces, conjugate flows and entropy production in this 
new language in sec |21 we go on to study the relaxation of a thermodynamic system 
in section IHl In particular we prove that if the shortest path principle is valid, then 
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the Universal Criterion of Evolution, written in a covariant form, is automatically 
satisfied. We also demonstrate that the term expressing the Universal Criterion of 
Evolution, written in a covariant form, has a local minimum at the geodesic line. 
Physically this means that a thermodynamic system evolves towards a steady-state 
with the least possible dissipation. We shall refer to this proposal, together with its 
corollary the shortest path principle, as the Minimum Dissipation Principle. 

In section 0] we obtain the thermodynamic field equations. For a — a ot j3 — (3 
processes we obtain the same equations found in refs but for general processes, 

we propose a new set of the field equations, which are conjectured to be also valid 
beyond the weak-field regime. These equations are notably simpler than the ones 
found in ref. 0]. In sectional we establish the equations that determine the steady- 
states and the critical values of the control parameters of a generic thermodynamic 
system at which a steady state becomes marginally stable. These equations are also 
conjectured to hold in the strong-field regime and therefore generalize the weak-field 
equations reported in refs Examples of calculations of the geometric quantities 

used in the new formulation of TFT can be found in the appendices. 

2 Forces, Flows and Entropy Production 

The central object of study in Thermodynamic Field Theory (TFT) is the ther- 
modynamic space, which is a smooth, path-connected, geodesically complete, real 
manifold M equipped with a Riemannian metric g. A Riemannian metric is a positive- 
definite quadratic form on a manifold's tangent spaces T^M, which varies smoothly 
from a point x of the manifold M to another. The metric g yields a smoothly-varying 
inner product 

{A, B) = A^g^,B\ A, Be T,M (2.1) 

where following the Einstein summation convention repeated indices are summed, as 
they will be in the remainder of this note. The inner product allows one to define the 
lengths of curves, angles and volumes. 

M contains one special point named x = 0, which corresponds to the unique 
equilibrium state. Note however that in practice the thermodynamical forces, which 
provide local coordinates on patches of M, will be not be independent but rather 
will satisfy a set of constraints that determines a submanifold N of M. When the 
submanifold N of physical configurations does not include equilibrium the system 
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relaxes do a different steady-state. 

Given two points a: G M and y G M, we may construct a shortest path 7(t) such 
that 

7 : [ti,tf] M, j{U) = X, j{tf) = y (2.2) 

The path 7 is automatically a geodesic. Clearly for some manifolds and choices of x 
and y there will be multiple paths 7 that minimize the length, thus it is crucial that 
observable quantities be independent of the choice of path. The metric g^y allows to 
define the invariant ajfine parameter r: 

c/r^ = gf.^dx^dx" (2.3) 

and our manifold M is "metrized" defining the length of an arc L by the formula 

L= [ \g^.x^x-y/^dT (2.4) 

The positive definiteness of the matrix g^^ ensures that rfr^ > and allows us to 
choose an affine parameter r that increases monotonically as the thermodynamic 
system evolves to a stable state. 

We introduce now the exponential map which is a map from the tangent space 
Tj-M of a Riemannian manifold M at the point x to another point exp^(V^) in M 
(see for example ref. [l^)- The point exp^{V) is the point ^vitf) on the geodesic 
7 that starts at 7y(ti) = x, where it is tangent to 7\/(ti) = V . tj is determined by 
the condition that the arc length along 7 from x to exp^.(\^) = 7y(t/) is the norm 
(V, V)^'^ of the vector V. This map is usually denoted by 

exp,(V) = 7y(t/) (2.5) 

We next introduce the generalized thermodynamic force vector Ux G T^M at each 
point X G M, which is defined along a shortest path 7[/^ by the relation 

exp,(t/)=7t/.(t/) = (2.6) 

Therefore, the exponential map Eq. ()2.6p associates to every tangent vector Ux G T^M 
the equilibrium point Ju^^itf) = in M. A thermodynamic system at x reaches the 
equilibrium state after traveling along the shortest path with velocity proportional 
to Ux at the point x. The norm | Ux \ indicates the geodesic distance from the point 
X to the equilibrium point measured along the shortest path 7^7^. The components U!^ 
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of Ux are interpreted as the generalized thermodynamic forces. The exponential map 
is a one-to-one correspondence between a ball in T^M to a neighborhood of x G M 
and so Eq. ()2.6|) defines a correspondence between the generalized thermodynamic 
forces and the points of manifold M. From now on we shall omit the suffix x being 
implicitly understood the dependence of vector U on point x. In Appendix A and 



[Appendix B| we will provide two examples of calculations of the map exp^{U). In the 
first example the thermodynamic space M is fiat space and so matrix tensor coincides 
with the Onsager matrix, while in the second M is the two-sphere. 

The thermodynamic flows J^, which are conjugate to the thermodynamic forces 
If^, are defined by 

= K^U" (2.7) 

where Xf^i, is an asymmetric tensor and, for brevity, we have suppressed the depen- 
dence of J, g and U on the point x and geodesic 7. Any 2-tensor may be decomposed 
into a symmetric and antisymmetric piece. The symmetric piece of A is identified 
with the metric tensor g and we shall name the skew-symmetric piece f^^^. Eq. ()2.7|1 
can then be re-written as 

J, = X.M" = {9,u + f,u)U' (2.8) 

is a one form J G T*M which, like [/, only vanishes at the equilibrium point 
X = although there may be points at which both U and J depend on a choice of 7. 
Eq. (j2.8|) allows one to generalize the inner product (A, B) of two vectors A and B 
in T^M from Eq. (j2IH) to 

{A,B)^{g,, + f,M'B'' (2.9) 



Finally we will define the entropy production a at the point x to be 

a = {U, U) = {g^, + f^.WU" = g^.M^U" > (2.10) 

where the inequality corresponds to the second law of thermodynamics. As the metric 
tensor (7^,^ is a positive definite matrix, this inequality is always satisfied. The entropy 
production a is an observable quantity and so the consistency of the formulation 
demands that it be independent of the arbitrary choice of 7. In fact a is just the 
square of the length of the shortest geodesic from x to 0, and so it is independent 
of the choice among geodesies of equal length and thus is globally-defined as was 
required. 
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In this note we will restrict our attention to homogeneous systems, therefore points 
X G M will correspond to configurations of the entire system. In particular the entropy 
fiux through the boundaries, will be constant. This implies that any variation in the 
entropy production a is equal to the variation of the time derivative of the total 
entropy. The generalization to inhomogeneous systems is straightforward and follows 
the strategy employed in Ref. (sf. 

3 Relaxation to a Steady-State 

In the TFT description of a thermodynamic system a homogeneous configura- 
tion corresponds to a point x in the thermodynamic space M. The corresponding 
thermodynamic forces, which are completely determined by the configuration a;, are 
assembled into the vector U defined in Eq. ()2.6|) . If the point x is not a steady-state 
(see section then the system will evolve. In this section we will consider the process 
known as relaxation in which the initial velocity vanishes and the system relaxes to a 
steady-state y. 

The TFT description of relaxation rests upon two conjectures 

a) The manifold (M,g) is a Riemannian manifold with metric g; 

b) The Minimum Dissipation Principle. 

The Minimum Dissipation Principle will be formulated in Subsec. where we 

will see that it implies that during the process of relaxation the configuration traces 
out a geodesic 7 in thermodynamic space. In Refs. the geodesic property was 

referred to as the shortest path principle. It allows one to write 

x = 7(0), y = ^{tf)=exp,iV) (3.1) 

where is a tangent vector to 7 at x whose norm is the geodesic distance from x to y. 
The system begins at x at time and reaches the steady-state y at the time tf. The 
process of relaxation begins as soon as the system is released, that is, immediately 
after time zero. In support of this conjecture, we will now argue that V automati- 
cally satisfies a manifestly covariant form of the Minimum Entropy Production, the 
Glansdorff-Prigogine Universal Criterion of Evolution and the Minimum Dissipation 
Principle. 
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3.1 The Minimum Entropy Production Theorem 

In 1945-1947, I. Prigogine proved the Minimum Entropy Production Theorem jlfl| . 
which concerns the relaxation of thermodynamic systems near equihbrium. This 
theorem states that: 

Minimum Entropy Production Theorem: Regardless of the type of processes 
considered, a thermodynamic system, near equilibrium, relaxes towards a steady-state 
y in such a way that the inequality 

^^(V,Vn^±iV,V)<0 (3.2) 

is satisfied during throughout the evolution and is only saturated at y. The Minimum 
Entropy Production Theorem has been proven by I. Prigogine, using a purely ther- 
modynamical approach. We shall now demonstrate that this theorem is automatic in 
the framework of TFT. 

Consider a geodesic 7(r) parametrized by r G [T"i,r/] equipped with a tangent 
vector field W normalized to have unit length 

W = ^l{r)=l*{^), {W,W) = 1 (3.3) 

Notice that r is not necessarily the time t, as the velocity is not necessarily unity 
throughout the evolution. The vector field W evaluated at any point x on the curve 
7 is a vector W{x) G Tr7 C T^M. We consider parameterizations of 7 such that 
W{x) 7^ 0. We then have the following identity 

Cwi,ir)){V,V) = ^{V.V^) (3.4) 

where Cw{-y{t)) denotes the Lie derivative with respect to the vector field W and, since 
we are near equilibrium, indices are raised and lowered using the Onsager matrix 
which is a metric for Euclidean space. 

There always exists a neighborhood U G M of y, which physically corresponds 
to the near-equilibrium approximation when y is near equilibrium, such that we may 
uniquely identify the coordinates x'^ of a point x with the vectors V{x) in a subset 
of a tangent fiber. One example of such a f/ is a ball whose radius is the injectivity 
radius, which in topologically trivial situations is roughly the inverse square root of 
the curvature, but any U such that each x and y in U are connected by a unique 
geodesic will suffice. To show that there is a correspondence between points and 
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vectors we need to show that given x we can determine V{x) and vice versa. We 
have aheady seen that we may uniquely define V{x) from the unique geodesic in U 
connecting x and y. 

To choose a name for the coordinate x given a choice of V{x) we must first confront 
the fact that V{x) is an element of the tangent space at x, T^M, and so to compare 
different V(x)'s we must somehow place them in the same space. This may be done 
canonically in U using the fact that in U there is a unique geodesic from x to the 
steady-state y. Thus we may define the vector Ay{x) G TyM by parallel transporting 
V{x) from X to y along 7. We may now use any choice of coordinates for y to write 
X on the same coordinate chart 

y-x = Ay{x) = V{x) (3.5) 

where we have used parallel transport to define a coordinate system for T^M which 
allows us to identify V{x) and A. Notice that each V{x) arises from a unique x because 
V{x) is parallel to 7 and the geodesic equation implies that parallel vectors remain 
parallel under parallel transport, and so A is also parallel to 7. Thus the uniqueness 
of the geodesic implies the uniqueness of A and so by Eq. ()3.5p the uniqueness of x. 
The sign in Eq. ()3.5|) is a result of the fact that V{x) points from x to the point y. In 
these coordinates the metric is the identity at y, the Christoffel symbols vanish at y, 
the geodesic curves in U that pass through y are straight lines in U and the entropy 
production is given by 

a=\x\^ + 0{x^) (3.6) 

Let us now evaluate the Lie derivative of {V, V) when the thermodynamic system 
relaxes from the point x to the stationary point y = 'y{tf). According to TFT this 
relaxation follows the unique geodesic in U that connects them, which is the straight 
line 

x = y + V{T), V{T) = V{x)f{T) (3.7) 
in the coordinates ()3.5p . We then find 

Cwi,ir)){V{T), V{t)) = 2{V{t), ^V{t)) = 2f{T)f{T){V{x), V{x)) < (3.8) 

where the last inequahty is a consequence of the fact that {V{x), V{x)) is positive by 
the positive definiteness of the Onsager tensor g, while /(r) is positive and /'(t) is 
negative because /(r) is the length of V which monotonically falls to zero during the 
relaxation. Only at the steady-state y is /(tj) = and so the inequality is saturated. 



8 



3.2 The Universal Criterion of Evolution 



The minimum entropy production theorem is generally not satisfied far from equilib- 
rium. However in 1964 P. Glansdorff and I. Prigogine demonstrated that a similar 
theorem continues to hold for any relaxation to a steady-state. We report the origi- 
nal version of the Glansdorff and Prigogine theorem [ll|: When the thermodynamic 
forces and conjugate flows are related by a generic asymmetric tensor, regardless of 
the type of processes, for time independent boundary conditions, a thermodynamic 
system, even in strong non- equilibrium conditions, relaxes to a steady-state in such a 
way that the following Universal Criterion of Evolution is satisfied 

Again Glansdorff and Prigogine have demonstrated this theorem using a purely ther- 
modynamical approach. In this subsection we shall prove that the Shortest Path 
Principle implies that a covariant form of the Universal Criterion of Evolution is 
automatically satisfied for relaxation processes in TFT. Consider a geodesic 7(t) 
parametrized by r G [tj, t/] equipped with a tangent vector field W 

W = -^i{t)=^*{—), {W,W) = l (3.10) 

Again we have parametrized 7 so that the velocity vector W always has a unit norm. 
As the velocity of the relaxation is an observable and is not necessarily equal to one, 
we must stress that r refers to a particular parametrization and not necessarily to 
time. However it is monotonically increasing with respect to time. As usual we also 
define a second vector field V on the curve 7 such that at any point x the vector V{x) 
satisfies 

7(ry)=exp,(\/(x)) (3.11) 

We will consider curves 7 which contain no closed loops, that is curves for which 
7(r) = 7(r') implies r = r' and so the condition ()3.11jl uniquely defines the vector 
field V. 

Now, our task is to evaluate the term in Eq. ()3.9j) when the thermodynamic system 
relaxes towards a steady-state along a geodesic. This term contains a time derivative 
of the vector V , which a priori is ill-defined as V at different times inhabits the 
tangent fibers over different points of M, which are not canonically isomorphic. The 
time derivative may be defined using the Lie derivative with respect to the vector W . 
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This leads to the Covariant Form of the Universal Criterion of Evolution 



(V, LwV) < (with {V, CwV) = only at the steady state) (3.12) 

From now on, the term (V, CwV) will be referred to as the Glansdorff-Prigogine 
quantity. Notice that V{'y{Tf)) = and so in particular the Glansdorff-Prigogine 
quantity relaxes to zero 

{V{^{Tf)), Cwi,ir,))V{^iTf)) = 0. (3.13) 

The quantity defined by the Lie derivative differs from the quantity in the mini- 
mum entropy production theorem, in fact the difference between these two terms is 
similar to the difference between the minimum entropy production term and the non- 
covariant time derivative term. The decomposition for the noncovariant derivative 
is 

d dV^ d 

= {V,^) + V^^[{g,. + f,.)V'^] (3.14) 

and for the Lie derivative it is 

^{V, V) = Cw{V, V) = {V, CwV) + V^{Cw{g + i)V), (3-15) 
dr 

We shall now prove the covariant form of the Universal Criterion of Evolution. 

Theorem: For all t < Tf we have the strict inequality {V,CwV) < 0. 

Proof: For every point x on the curve 7 the vectors V{x) and W{x) are elements of 
the one-dimensional vector space T^.^, therefore they are parallel and there exists a 
function fipfir)) such that 

V{l{T)) = f{^{t))W{^{T)). (3.16) 
We may now evaluate the Lie derivative 

Cwi,ir))V{^{T)) = [iy(7(r)),\/(7(r))] = [iy(7(r)),/(7(r))iy(7(r))] 

d d dfh(T)) 

= /(7(r))7*(^)] = -^^W^(7(r)). (3.17) 

Combining Eqs. ()3.16|) and ()3.17|) we may evaluate {V, CwV) as a function of / 

{V,CwV) = ifW^^W) = f^{W,W). (3.18) 
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The fact that g is Riemannian and so positive definite imphes that {W, W) > 0. We 
have seen that V^(7(t/)) = and we have imposed that W therefore /(tj) = 0. 
In addition the fact that W ^ imphes that the distance from 7(r) to 7(tj) has no 
critical points and so / has no critical points. Thus / is strictly monotonic in r and 
at the maximal value t = Tf it vanishes. As a result / and df/dr have opposite signs 

< (3.19) 

and so the right hand side of Eq. ()3.18p is the product of a positive and a negative 
term. Therefore (V, CwV) is negative as desired. 

At this point one may object that the Glansdorff-Prigogine quantity in Eq. (j3.12p 
is not actually the quantity that appears in ()3.9|1 because t and r are not equal and 
so W is not actually the velocity of the system. This implies that the covariant form 
of the Universal Criterion of Evolution is not in general equivalent to the original 
form. However after a little algebra one can see that they are equivalent whenever 
the deceleration \a\, velocity v and affine length / of the system satisfy the following 
inequality 

\a\ < (3.20) 

/ — r 

at every affine time r. The acceleration is unconstrained, and so this is only a 
constraint on the dissipation of the system as it approaches the steady-state and 
slows. However near the steady-state the denominator of the right hand side of 
Eq. (j3.2Up shrinks to zero and so even near the end of the relaxation the constraint 
on the rate of dissipation appears to be quite weak. This constraint will be compared 
with experiment in a subsequent publication. 



3.3 The Minimum Dissipation Principle 

In refs ^-0 we have introduced the Shortest Path Principle. In Subsec 13. 2[ we 
have demonstrated that if the shortest path principle is valid then the Universal 
Criterion of Evolution, written in a covariant form, is automatically satisfied. In this 
subsection we will show that the expression (V, CwV) has a local minimum at the 
geodesic line. We refer to this theorem, together with the Shortest Path Principle, as 
the Minimum Dissipation Principle. The Minimum Dissipation Principle corresponds 
to the following physical conjecture: 



11 



For any type of physical process, even in strong n on- equilibrium conditions, a thermo- 
dynamic system subjected to time-independent boundary conditions relaxes towards a 
steady-state following the path with the least possible dissipation. 



Let us consider an arbitrary path 7(t) parameterized by t G [0,tf] equipped with 
a tangent vector field 

This time we fix the velocity of the trajectory, this is necessary in any minimization 
of dissipation because otherwise one could always trivially speed up the dissipation 
by accelerating the process. In particular we will parametrize 7(t) by imposing that 
is a unit vector 

{W,W) = 1 (3.22) 
At every time t there exists at least one geodesic Q;t(r) such that 

«i(0)=7(t), at{Tj)=^itf) (3.23) 

There also exists a vector Vt in the tangent space T^(^t)Oit C T^(^t)M such that 

7(t/)=exp,(,)(\4) (3.24) 

Eq. ()3.24|1 defines the vector Vt at every point 7(t) on the curve 7 and so it defines a 
vector field V{'-f{t)) on 7 with vectors in the tangent space TM. Thus we may define 
the Lie derivative Cw(yit))V{l(t)) which will be another vector field on 7. 

We will show that, using the above parametrization 

{V, CwV) = -\V\cos{e) (3.25) 

where 6 is the angle between a and 7. This quantity is minimized when = 0, 
corresponding to the case in which a and 7 are parallel. Therefore the Glansdorff- 
Prigogine quantity is minimized everywhere iff a and 7 coincide everywhere, in which 
case 7 is a geodesic. 

To prove this theorem we shall perform our calculations in Fermi coordinates 
jisl . In particular, we shall work in the Fermi coordinate system such that the first 
derivatives of the components of the metric along the geodesic at vanish 

^ g^uiat) = (3.26) 
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Figure 1: As one moves along the curve 7 the minimal distance to 7(t/) shrinks 
at a speed of cos{9). This speed, which is proportional to the negative quantity 
{V,CwV), is locally minimized when 6 = corresponding to the case in which 7 
follows a geodesic. 



Therefore in Fermi coordinates, a distance e from the geodesic, the components of 
the metric will differ by a correction of order e^: 

- + + + Oie') (3,27, 

or 

Ag,. = gA<^t+e) - g^o^t) = ^^e^^'^ + O(e') (3-28) 

In Fermi coordinates, we may interpret the Glansdorff-Prigogine quantity ()3.25j) using 
the fact that (V^, V) is the length squared of the geodesic at and its Lie derivative 
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with respect to W is 

Cw{V, V) = 2{V, CwV) + V^V\Lwz)^u (3.29) 

Notice that the term Cw^-, ^) (^oes not correspond to the Lie derivative of the entropy 
production, which as seen in the section |21 is defined to be Cw{U, U). 

We shall now show that the second term on the right-hand side is of order e while 
the other two terms are of order one 

V^ViCw^)^, ~ 0(e), Cw{V, V) ~ 2{V, CwV) ~ 0(1) (3.30) 

Up to now, we have defined the vector W only along the curve 7. However, Eq. (|3.29p 
requires a definition of the derivative of W in directions that are not parallel to the 
curve. Such derivatives have not yet been defined, but they may be defined if we 
extend the definition of to a neighborhood U of the point at(0). We define W at 
a point at(0) + e'^ in this neighborhood by parallel transporting Wjati^)) along the 
unique geodesic 5x^ = 'm.U that connects aj(0) and 0:4(0) + e'*' |lj| (see Fig. Q): 



W(Y(t>fe) 




An Arbitrary Path Y 



Figure 2: Definition define the vector W in the neighborhood of the point 7(t) by 
parallel transporting the vector W{''){t)) along an edge 5x^ = e^. 

Wtix) = Wi; + T^,^{x)e'W,^ + 0{e') (3.31) 
14 



Eq. (j3.3H) allows us to generalize the definition of the quantity (V, C-w'^) to arbitrary 
curves. The last term of Eq. ()3.29p is the derivative of the metric along the direction 
parallel to a^. It is easy to check that this term is of order e in our coordinate system: 



V^V\L^Z)p.u = y^V^W^-^ + '^V.V^-^ (3.32) 



or 

V^V^iCw^),, = 2\/^\/^iy^rC, + 2V^V''^e''WS + 0{^) (3.33) 
where we have used the identity 

^ = n,9.u + n^g., (3.34) 

However in Fermi coordinates, a distance e from the curve the components of the 
affine connection will differ by a correction of order e: 



rj.(a.+.) = TlMt) + -^^^^ + Oie') = -^e- + ^(e^) (3.35) 



Observe that from Eq. (j3.35j) . we also have that 

^ = lim ^ ' ~ 0(1) (3.36) 

Therefore, the second term on the right-hand side of Eq. ()3.28p is of order e. As 
a check, we note that Eq. ()3.3(i|l reveals that also the Riemannian curvature tensor 
R'^Xk is of order 0(1), as it must be since our choice of coordinates cannot affect the 
curvature. 

Now, we shall prove that the derivative of the length squared of at along the 
curve 7 is of order 0(1). For this, we have to compare the lengths of the geodesies 
at and at+e- This leads us to the notion of Jacobi fields (see, for example, ref. 
Let us consider a single infinite family of geodesies in the manifold M. Let r be a 
parameter varying along each geodesic of the family, and let t a parameter constant 
along each geodesic of the family, but varying as we pass from a geodesic to another. 
A Jacobi field is a vector field on a geodesic that describes infinitesimal deformations 
that interpolate between different geodesies in a one-parameter family. For example, 
the family at yields the Jacobi field Jt 

JtVr) = —Qf- (3-37) 

on each geodesic (r). 
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Consider now two nearby geodesies in the a family, at and at+e- We call the 
points A on at and N on at+e corresponding points if they have equal values of r (see 
ig. El). Let = esin^Jf be the infinitesimal Jacobi vector (see Fig. H} such that 



g^^a^tK = (3.38) 

This equation tells us that the deviation is perpendicular to a^. The vector Xt 
satisfies the Jacobi equation 3| 

;^Ar(r) + K,^^^XKr)^y- = (3.39) 

where D denotes the covariant derivative. Thus both terms in the Jacobi equation 
()3.39|1 are linear in e. This implies that, in the limit e — 0, the second derivative 
of our Jacobi field goes to zero. The relative angle at the point A' (see Fig. 
between at and at+e, measured using parallel transport along the flow of the Jacobi 
field, can be easily evaluated taking into account the relation 0] 

r f da^ 
{/\a^t)AY,)^, = j j Y,R^^,^a'',^\^dtdT (3.40) 

where 1^ is a vector chosen arbitrarily at A'. Eq. ()3.40|) shows then that (p is of order 
e. 

Our task is to compute the hmit of quantity 

Ito = 2s^V' lim ^ + V"V^ lim ^ (3.41) 

Using the relation ()3.28|) . the last term on the right-hand side vanishes. As seen in 
Fig. El to leading order in e, the difference between the vector at the geodesic at 
and its image on the curve at+e under expat{Gsm{9)Jt) is of order e^. In conclusion, 
the limit of the difference between the squared length of the geodesic at+e and the 
portion of the geodesic at such that at+e = exp^(()(esm(6') J^), vanishes. 

However, the curve at+e is not quite the same as expat{esin{9)Jt), because we have 
parameterized at+e such that (see Fig. IH) 

at+e{r = 0) = 7(t + e) ^ exp^^,^{estn{e)Jt{0)). (3.42) 

That is to say that the Jacobi field is not parallel to 7, so it does not take the beginning 
of the curve at to the beginning of the curve at+e- Instead exp of the Jacobi field 
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Figure 3: When the relative angle (j) at the point A' between at and at+t is of order 
e, the difference between length A'B' and length AB is of order e^. Indeed, we have 
A^B' - AB ~ \(I)CB' ~ 0{e^). 

esin(^) Jt(0) takes the beginning of the curve at to a point — ecos(^) on the curve at+e- 
Up to corrections of order the curves at and at+t are the same length, and so the 
distance from at+eiS^) to 7(t/) is that of at minus the correction ecos(6') that we have 
just calculated. The distance squared then differs by |V^|ecos(6'). By definition the 
Lie derivative is this quantity divided by e, which yields the claim (j3.25j] . 

4 The Thermodynamic Field Equations 

As seen in Eq. ()2.10|) . the second principle of thermodynamics is satisfied if the 
metric g is positive definite. However the postulates that we have presented so far 
are insufficient to calculate g. In this section we will add a third postulate to TFT 
stating that the metric is the solution to a set of field equations which are in turn 
equal to the variations of an action functional. In short we include the postulate: 

c) The Principle of Least Action. 

In particular we will assert that the action is invariant with respect to general coordi- 
nate transformations and we will not include a cosmological constant term, as there 
has been no experimental evidence for the presence of such a term in TFT. 
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We will also place the following three restrictions on our physical systems 

1) Internal fluctuations are neglected; 

2) Boundary conditions are time-independent; 

3) The system is not subjected to external perturbations. 

If these three conditions are simultaneously satisfied then the equations of motion 
will be source-free. The internal fluctuations can be analyzed using the Landau- 
Lifshitz theory and systems submitted to time-dependent boundary conditions have 
been examined by York in ref. Examples of systems subjected to external 

perturbations, those most commonly encountered in practice, can be found in ref. 



4.1 a-a and f3-/3 Processes 

In a-a and jS-jS processes the matrix X^^, that relates the thermodynamic forces and 
flows ()2.7|1 is symmetric and so is equal to the metric g. We therefore expect the action 
to be constructed entirely from g. After the cosmological constant, the invariant 
constructed using the least number of derivatives of g is the Ricci scalar R = g'^^R^^ 
where R^j,^ is the Ricci tensor [3| 

p _ pA _ pA I pr? pA _ pT? pA (A^'\ 

Of course one could add higher powers of R and its derivatives to the action and 
maintain general covariance, but so far Eq. ()4.H1 has been sufficient to reproduce the 
experimental data and numerical simulations. If such terms are indeed present in the 
action then, by dimensional analysis, they will be multiplied by characteristic scales of 
positive spatial dimension and they will be negligible at smaller distance scales. Thus 
it may well be that above some threshold length scale, corresponding to very strong 
fields, higher order terms in R will begin to play a role. For example in gravitational 
theories such terms may play a role deep inside of black holes and immediately after 
the big bang, but are irrelevant at the scales of all observed phenomena. By contrast 
in general relativity the nonlinearity of the R term considered here already provides 
a significant correction to the clocks of the GPS satellites and its effect on Mercury's 
orbit was observed a hundred years ago. 
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We therefore postulate that the action / takes the simple form 

1 = j ^R^^g^^dn (4.2) 

where dVt is a volume element and the domain of integration is the entire manifold 
M. Here g denotes the determinant of the metric g^j_p. The principle of least action 
implies that the variation of I with respect to the metric vanishes, which yields the 
familiar thermodynamic field equations: 

R^.u - ^Rgf^u = (4.3) 
For n 7^ 2 Einstein's equation (j4.3p reduces to the Ricci flatness condition 

R>.. = (4.4) 



4.2 The General Case 

Let us now consider the general case in which A^,y is not necessarily symmetric 

= Vf^' = i9f.u + MU" (4.5) 

Any completely skew-symmetric tensor of type (0, 2) defines a 2-form /. Given a skew- 
symmetric tensor / of type (2, 0) with components f^'^, we can define an (n — 2)-form 
*/ called the Hodge dual of /. In terms of components, the relations between */ and 
/is ^ 

where £ij.ujij2...j„-2 is the Levi-Civita symbol: 



_ i,j...n 

-ij...n 



+1 if ij .. .n is an even permutation of 1, 2 . . . , n 

— 1 if zj ... n is an odd permutation of 1, 2 . . . , n (4.7) 

otherwise 



The inverse of Eq. fl4.6|) reads 

To write the kinetic term for the /-field we will need to first define its antisym- 
metric 3-tensor field strength H 

H^up = d^f^p + d^fp^ + dpfp,^ (4.9) 
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or simply 

H = df (4.10) 

where d is the exterior derivative. The square of the exterior derivative, d"^ contains 
the contraction of two ordinary derivatives with the e tensor and so it vanishes by 
antisymmetry. Therefore H is closed 

dH = d'^f = (4.11) 

In principle the balance equation for H could also contain a source term J' 

dH = J (4.12) 

where JT" is a closed 4-form, dj' = 0. 

Given a third-rank skew-symmetric tensor H we can define an [n — 3)-form *H, 
the Hodge dual of H, with components 

*-^iii2-Jn-3 = '^^9^^'^^t^^^hi2-i-a-zH^^'^ (4-13) 

As we shall see shortly, *H also satisfies a balance equation with a source term, which 
we shall denote *J. This source term may be nonzero only if the three restrictions 
on our physical system are relaxed. In other words internal fluctuations, external 
perturbations and time-dependent boundary conditions can serve as *J sources for 
*H. However in all cases examined so far no potential J sources have been found for 
H. 



We are therefore led to the following postulate [3|,|^: 

d) The 3-form H has no sources or, equivalently, the 3-form H is closed. 

Notice that, from the mathematical point of view, to impose J7 = it is sufficient to 
require that the tensor field / does not possess any singularities and so Eq. ()4.10|1 is 
well-defined everywhere. 

We will now present an invariant action / for the thermodynamic system. As in 
the previous subsection, we will restrict our attention to terms that contain the least 
number of derivatives, as higher derivate terms will be suppressed by a characteristic 
length scale and so may only be relevant in the presence of very strong fields. Again 
the metric may enter via an Einstein-Hilbert term equal to the Ricci scalar R. The 
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antisymmetric tensor / may enter via a 2-derivative term H A *H generalizing the 
Maxwell action and also via a 0-derivative mass term xf^*f where x is a dimensionful 
parameter. Here A is the wedge product of differential forms and * is the Hodge star 
defined in Eqs. ()4.6|1 and ()4.13j) . A cosmological constant term and more complicated 
functions of the mass term are not necessary to explain the data and so will not be 
included. In components the proposed Lagrangian density is then 

Ck^n = ^H,,,H^^^ + ^f,,r (4.14) 

where indices are raised using g. The mass parameter x is an inverse of the charac- 
teristic length squared of the thermodynamic system. 

The source J^^ can be included to the theory by adding an interaction term to 
the Lagrangian density 

C^nt = Sn-lJ,uF' (4.15) 

where Sn-i = 2'K^/'^V{n/2) is the area of a unit [n — l)-sphere in n space dimensions 
{V{n/2) is the Gamma function). Notice that if we do not add the source J, we could 
add a new term coupling the source J to the dual field /('^"'^') where ^/('^""O = 
However, / and are not mutually local, since at a point cannot be 

constructed from / at that point together with a finite number of derivatives and so 
no local theory can simultaneously contain both J and J. In line with our postulate, 
we impose JT" = and so we only consider the source term in Eq. ()4.15|1 . 

Combining ()4.14|1 . ()4.15j) and the Einstein-Hilbert action we arrive at the total 
action 

1 = J^C = J i^R +^HA^H+ |/ A + J A */) (4.16) 

The relative normalization of the first two terms in the action (|4.16p may be 
modified by rescaling the metric, here we have chosen the convention known as the 
string frame in string theory |2]] |. 

*J will be treated as an external source, and so there are two variables whose 
variations yield the field equations. The variation with respect to the metric yields a 
modified Einstein's equation 

Rfiiy- 291^^1^^ = ~ I iH^XpH^^--gf_,uHxK.pH^'^'')-x(^fiiKf,y'^-^9iiiyf\K.f^'^^ H — "^^^^^.u-hnf^'^ 

(4.17) 
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while the variation with respect to the /-field yields 



d* H + X* f = ~Sn-i * J 



(4.18) 



We have seen that the antisymmetric tensor / modifies the field equations for the 
metric. One may wonder if it also affects the shortest path principle. In general a 
p-tensor potential can be covariantly coupled to a p-dimensional J-source and to a 
{d — p — 2) -dimensional ^/-source. For example in 4-dimensional Maxwell theory the 
1-tensor potential couples to electric and magnetic particles, while in 5 dimensions 
the electric sources are still particles but the magnetic monopoles are strings. In our 
case p = 2, so the J-source is a string and the jT-source is {d — 4)-dimensional. In 
particular, the trajectories in the manifold of thermodynamic configurations are, being 
paths, 1-dimensional and so they cannot contain the J-source. Except in 5-dimensions 
they also cannot contain the jT-source. Thus, there is no covariant coupling of the 
trajectories to the /-field, as there would have been for an ordinary vector potential. 

Concretely, this means that the action for a relaxing configuration is proportional 
to the proper length of its trajectory 



just as in the case with no /-field. The proper length of a path, by definition, is 
extremized by a geodesic. Thus the solutions of the classical equations of motion 
for the evolution of a state are still geodesies in thermodynamic space, even when 
we include the / field. This is not to say that the solutions with an / field are just 
the solutions to the system without the /-field, because the new Einstein equation 
()4.17|1 no longer admits Ricci-fiat metrics, and so the Riemannian manifold M itself 
is changed by the inclusion of /. 

5 Steady-States and Stability Criteria 

In the thermodynamic theory of irreversible processes, a steady-state of order k, 
where k is less than or equal to the dimension n of our thermodynamic space M, is 
a thermodynamic point x with coordinates 




(4.19) 



U\ U\..U' 



kept constant 




a 
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From now on, Latin indices will run from k + 1 to n and Greek indices will run from 
1 to n. Thus fj5.1|) may be written more compactly as 



J>.dU^ \st.state= (5.2) 

where symbol d denotes the total differential on the submanifold N of solutions of 
the constraint equations. 

Steady-states may be stable or unstable. According to the Layapunov stability 
theory the stability of a state is determined as follows 

if |r/M < then Un is unstable for t > tn 

dt ' " ~ 

if 'ji^'^S Ij/m > then Uq is asymptotically stable for t > to (5.3) 



dt 



where S denotes the total entropy of the thermodynamic system. In practice we can 
modify the value of by varying a series of control parameters A. When these 
parameters reach a critical value Ac, the sign of the inequality ()5.3|1 will change, and 
a steady-state will lose its stability. A steady-state in a configuration with critical 
control parameters is referred to as a state with marginal stability: 

j^S^S 1^. (Ac) = for t > to (5.4) 

The total entropy S is related to the entropy production a considered in this note 
via the entropy balance equation 

^S = a + Es (5.5) 
dt 

where S5 denotes the entropy flux through the boundaries. Eq. (j5.5p may be varied 
twice to obtain 

^S^S = 6^a + 6^J:s (5.6) 

where the variation of the S5 term vanishes because we consider homogeneous config- 
urations in which S5 is constant. In ref. ^|, the authors considered inhomogeneous 
thermodynamic systems subject to boundary conditions in which all fluxes vanish at 
infinity. We are instead considering finite homogeneous systems for which the en- 
tropy flux through the boundaries, is constant. Notice that in the particular case 
S5 = the only steady-state is equilibrium. In all other cases our systems do not 
satisfy their boundary conditions, as they are considering all of spacetime whereas we 
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consider only a region in which the system is homogeneous. As a result simplifications 
that occur in ref. ll| which rely upon their choice of boundary conditions do not 
apply here, but they will not be needed. Thus in our case the stability condition is 



if (5^0- |r/M < 

if 6^ a If/M > 



then Ug is unstable for t > to 



then Ug is asymptotically stable for t > tg 



(5.7) 



These conditions may be interpreted geometrically as follows. Recall that the 
point X is in an j dimensional submanifold C M of configurations that satisfy a 
set of physical constraints that relate the thermodynamic forces. Then x is a steady 
state of order k if j = n — k, and if the derivatives of a with respect to all tangent 
vectors to vanish. Diagonalizing the matrix of second derivatives we may also define 
the stable directions to be the positive eigenspace, the unstable directions to be the 
negative eigenspace, and the marginally stable directions to be the null eigenspace. 

The steady-state condition ()5.2|1 is not canonically defined as the f/'s at various 
points X live in different spaces. However using the geometric interpretation of steady- 
states as critical points of S one arrives at a unique covariant definition for steady- 
states, and thus a definition that we conjecture will hold beyond the weak-field regime. 
The steady states are critical points of a which means that da vanishes where d is the 
exterior derivative on the submanifold A^. As a is a scalar, the ordinary and covariant 
derivatives of a are equal and so 



where the / term vanishes because / is antisymmetric and UU is symmetric and Dg 
vanishes because the covariant derivative of the metric always vanishes. Thus the 
covariant form of the steady-state condition is just 



The stability condition (|5.7|) is already covariant, as is a scalar and so ordinary 
derivatives on S are canonically defined. But for completeness we evaluate the double 
variation 



= da = Da = DiU'^Wg,,,) = 2g^M''DU^ 




g^M'DU^ = 0. 



(5.9) 



(5.10) 



dx"-dx^ 
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where 

Q^,=29,Mi'aU:',+2f,.-,aU;^U;t+2ig,^+f^^^^^ (5.11) 
We have again used the fact that g is covariantly constant. 

There are marginally stable states when the determinant of the symmetric matrix 
alfe vanishes 







\=\c 

(5.12) 

We expect that, even beyond the weak-field approximation, Eqs ()5.9p and ()5.12|) 
will respectively determine the steady-states Uq and the critical values Ac of control 
parameter A. 



6 Conclusions 



We have presented a manifestly covariant and coordinate-independent formulation 
of the Thermodynamic Field Theory (TFT). This allows the theory to be extended, for 
example, to thermodynamic spaces M that cannot be covered by a single coordinate 
patch corresponding to theories for which different observable quantities are suitable 
in different regimes. Making use of the Shortest Path Principle, we have demonstrated 
the validity of the Universal Criterion of Evolution, expressed in a covariant form, 
and we have shown that this term has a local minimum at the geodesic line. We 
referred to this theorem, together with the Shortest Path Principle, as the Minimum 
Dissipation Principle. Physically, the Minimum Dissipation Principle affirms that, 
for time-independent boundary conditions, a thermodynamic system evolves towards 
a steady-state with the least possible dissipation. 

Of course the test of any proposed formulation is its agreement with experiment. 
One physical system for which a large amount of data is already available, and a 
theoretical underpinning is lacking, is the relaxation of a magnetically confined plasma 
towards a steady-state in the nonlinear classical and Pfirsch-Schliiter regimes or in the 
nonlinear Banana and Plateau regimes. The metrics in these regimes can be found 
in 0. 

We have also derived the thermodynamic field equations. For the a — a oi (3 — (3 
processes, these equations reduce to the ones already found in refs 0-01 • The valid- 
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ity of these equations, in the weak-field approximation, have been tested analyzing 
several thermodynamic systems. So far, we have found that numerical simulations 
and experiments are in agreement with the theoretical predictions of the TFT. When 
the skew-symmetric tensor /^i, does not vanish, we have proposed a new set of field 
equations, which are conjectured to hold even beyond the weak-field regime. The 
strong field equations are, however, notably simpler than the ones obtained in ref. 

The equations that determine the steady-states and the critical values of the 
control parameters of a generic thermodynamic system have also been presented. 
These equations are also conjectured to hold beyond the weak-field regime and so 
they generalize the weak-field equations reported in refs 0-j4|. We note that this 
extension of linear thermodynamics is quite different from the Prigogine-Glansdorff 
extension, which is based on hydro- (or plasma-) dynamical stability theory and 
bifurcation analysis 0|. The relation between the two approaches should be made 
explicit. This will be the subject of a further publication. 

The current approach describes the trajectory followed by a relaxing system, but 
does not give the time dependence of the relaxation process. It may be possible to 
include the time direction in the thermodynamic space with a metric such as gu = — o". 
The geodesic evolution in spacetime would then give a firm prediction for the time 
dependence of the spatial trajectory. In addition the steady-state condition and the 
Layapunov stability conditions may be derived from the constraint that particles 
follow spacetime geodesies. Of course in order for the system to stop at a steady 
state somehow dissipation will need to be added to the system, once this is done it 
may be possible to apply TFT to processes other than relaxation. 



Appendix A Example: Flat Space 

Let M = vr be a plane. We will compute the map exp : Ttt i-^ tt. In this case the 
Riemannian manifold (M, g) coincides with (M, L) where L is the Onsager matrix. 
There is a well-defined exponential map from a ball in T^tt to a neighborhood Ux 
of X G tt: 

C T,.7r = C7i = (A.l) 

When the metric is fiat, the geodesies 71/ (r) are straight lines 

= tA^" + (A.2) 
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where is a constant vector and A G T^ir, 'jA{ti) = x and 'jA{tf) = Xf. Therefore 
the exp map on x is simply 

exp^TA = tA + x (A.3) 

that is, the exponential map exp^ V (with p E n and V G T^tt) simply performs a 
translation: 

expp V =p + V (A.4) 

Let us now consider the case U E d TqM = B? where G vr denotes the origin, 
corresponding to the equilibrium state. In this case U points from x to the origin 
and, using the canonical identification of the plane vr and its tangent fibers, which 
are planes of the same dimension, we find the relation x = —U (see Eq. (13. 5j) ). If the 
coordinates of x are x = (a, b), the coordinates of U are then U = (—a, —b) and, from 
Eq. ()A.4|) . we have 

ex.ppU = p + U = {a,b) + {-a,-b) =0 (A.5) 

This implies that the exponential map sends the tangent vector U at p to the equi- 
librium state, which is the definition of U in Eq. ()2.(ij) . 

The entropy production of a system at the point x is just its length squared 

a = = I = (A.6) 



Appendix B Example: The Two-Sphere 

Let M be the unit two-sphere 5*^ with the usual rotation-invariant metric 

ds^ = de"^ + sm'^ed(f)'^. (B.l) 

We will use both spherical and Cartesian coordinate systems, which are related by 

x = sin6' COS0, y = sm9 sin0, z = cos9. (B.2) 

The geodesies on the 2-sphere are the great circles, which are the equators corre- 
sponding to various choices of north pole. We will now fix the north pole to be = 0, 
and so every geodesic 7 will intersect the equator = 7r/2 at two points = —6 and 
(J) = IT — 6. Note that the equator itself is the geodesic that intersects itself an infinite 
number of times, to apply the following derivation to the equator itself it would be 
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necessarily to deform it by a small rotation and then take that rotation to zero at 
the end. However the final result will make sense for every geodesic, including the 
equator. 

The map exp{v) is easy to calculate in the case in which 7 is the equator, it just 
increases (p by the norm v of the vector v. To evaluate exp on an arbitrary point 
{60, (po) with an arbitrary tangent vector v, which yields an arbitrary geodesic 7, it 
will be simplest to rotate the coordinate system so that 7 becomes the equator, then 
exp may be evaluated and the coordinate system may be rotated back to yield the 
final answer. 

The rotation of the coordinate system will proceed in two steps. First we use 6 to 
define a rotated set of coordinates 



a; = sin^ cos((/) + 5), y = sm9 sin(</) + (5), z = z = cos9. 



(B.3) 



Then we rotate about the y-axis until 7 is the equator, which we parametrize with 
the azimuthal variable r 
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sin a 


cos a 


/ 




(B.4) 

sin a cos a y y J 
In coordinates the rotation ()B.4|) is 
sin 6 cos((/) + 5) = cos r, sin 6' sin(</) + 5) = cos a sinr, cos = sin a sinr. (B.5) 



The tangent vector v will be decomposed in polar coordinates {v,(3) where v is 
the norm of v and (3 is the angle between v and the azimuthal direction 0. In this 
section we will calculate exp[0o,(j}o)V. First we will find expressions satisfied by the 
unknown angles r and a in terms of the given data 6*0, (j)o,v,j3. These will allow us 
to calculate the z component z' of the function exp, which is particularly simple as 
it is independent of 6. Then we will rotate the coordinates to construct the other 
components x' and y'. 

w is a tangent vector to the geodesic 7 and so the rotation ()B.4|) relates v to the 
equator's tangent vector d/dr. For example the last expression in ()B.5|) yields 



d d . 

— cos U = — sm a sm r 
OT or 



(B.6) 
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where the derivative is evaluated at (6*0,00)- The left hand side may be evaluated 
using the chain rule, where we note that the derivative of 9 with respect to r along 
the geodesic is just the theta component sin/? of the unit tangent vector v/v. The 
right hand side is evaluated by noting that a is the angle by which the entire geodesic 
7 is rotated, and so it is independent of r. Thus ()B.6|) reduces to 

sin/3 sin6'o = sinQ; cost. (B.7) 

The dependence on the unknown angle a may be removed using the last term in 
Eq. ()B.4j) . We then divide both sides of the resulting equation by cos^^o to obtain 

cotr = sin/5 tan6'o. (B-8) 

Thus we have determined the unknown angle r in terms of the initial point and 
tangent vector. While in principle we could also find a, it will suffice to recall the 
last term from ()B.5|) 

sin a = cos ^0 cscr. (B.9) 

In our new coordinates the exp operation corresponds to a simple shift of the 
equator's longitudinal coordinate r 

T^T + V. (B.IO) 

Finally we need to reexpress this shift in terms of the old coordinates. In fact it will 
suffice to find the corresponding coordinate z' . Using the rotation ()B.4jl at the new 
coordinate t + v 
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we find 

z' = sma sin(r + f). (B.12) 

This may be expanded using ()B.9|) and the trigonometric identity for the sine of a 
sum to yield 

2;' = cos 6*0 cscr(sinr cosf + cosr sint>) = cos6'o(cosf + cot r sint>). (B.13) 
Finally we may use ()B.8|) to obtain 

z' = cosOq cos f + sin/? sin 9q sinf. (B.14) 
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To find the corresponding formula for x' we may rotate tlie coordinates by 7r/2, 
yielding a new z coordinate, 5, equal to the old x coordinate 

5 = cos ^ = sin 6' cos0 = a;. (B.15) 

This rotation is an isometry so v is unchanged, however the angle /3 is effected. The 
rotated value of /3 may be found similarly to ()B.6|) 

d ~ ~ ~ a 

— cos ^ = sin/? sin0 = — — (sin6' cos </>) = sin /3 cos 5 cos(/) — cos/3 sin</). (B.16) 

OT OT 

In the last term of the last expression of ()B.16|) we have used 

cos/3 = ^/^—, ^<^<^ = sm^. (B.17) 

We may now use ()B.15|) and ()B.6|) to rotate our equation (jB.14|) for the z compo- 
nent of the exp function to obtain the x component 

x' = sin^o cos</)o cos f + (sin /5 cos 6*0 cos </)o — cos /? sin^o) sinw . (B.18) 
And at last we may rotate to obtain the %j component 

?/' = sin6'o sin 00 cosf + (sin/3 sin 6^0 cos 00 + cos /5 cos0o)sini;. (B.19) 
Putting together Eqs. ()B.14IB.18IB.19j) we obtain 

Taking the north pole to be equilibrium, the thermodynamic force at any point 
is calculated using the geodesic that runs straight north from any point to the north 
pole 

f/ = (t, = ^,/3 = 0). (B.21) 

This geodesic is uniquely determined for all points except for the north pole 9 = n, 
for which there are an infinite number of northward geodesies and the coordinate (3 
is ill-defined. The entropy production is then simply 

a = {U, U) = e'^ (B.22) 

Note that at the south pole, where the geodesic is not unique, a is independent of 
the geodesic chosen and so continues to be well defined. 



30 



Acknowledgments 

We acknowledge our indebtedness to Prof. CM. Becchi of the University of Genova 
- Physics Department and to Prof. E. Tirapegui of the University of Chile - Physics 
Department, for reading this manuscript and making helpful suggestions. 

One of us (GS) would like to express his sincere gratitude to his hierarchy at the 
European Commission: Prof. A. Mitsos, Dr P. Fernandez Ruiz and Dr E. Rille. He 
is also very grateful to Dr U. Finzi and Dr M. Cosyns, of the European Commission, 
for their continuing encouragement. 

GS would like to thank the members of the EURATOM Belgian State Fusion Asso- 
ciation, in particular Dr D. Carati, Dr B. Knaepen and Dr B. Weyssow. 

The work of JE is partially supported by IISN - Belgium (convention 4.4505.86), by 
the "Interuniversity Attraction Poles Programme - Belgian Science Policy" and by 
the European Commission RTN program HPRN-CT-00131, in which he is associated 
to K. U. Leuven. 

References 

[1] Sonnino G. // Nuovo Cimento, 115 B, (2000), 1057. 

[2] Sonnino G. 2001 Thermodynamic Field Theory (An Approach to Thermody- 
namics of Irreversible Processes), proceedings of the 9th International Work- 
shop on Instabilities and Nonequilibrium Structures, Vina del Mar (Chile). 
Kluwer Academic Publishers, Nonlinear Phenomena and Complex Systems 
Editions, 9, (2004) 291. 

[3] Sonnino G. 2002 A Field Theory Approach to Thermodynamics of Irreversible 
Processes, (These d'Habilitation a Diriger des Recherches - H.D.R.) - Institut 
Non Lineaire de Nice (France). 

[4] Sonnino G. // Nuovo Cimento, 118 B, (2003) 1175. 

[5] Sonnino G. Int. J. of Quantum Chemistry 98 (2004) 191. 

[6] Sonnino G. 2005 A Field Theory Approach to Transport Processes in Confined 
Plasmas: the Nonlinear Classical and Pfirsch-Schliiter Regimes and Transport 



31 



Equations in the Nonlinear Banana and Plateau Regimes (submitted to pub- 
lication in the review I.A.E.A. Nuclear Fusion). 

[7] Peeters Ph. and Sonnino G. // Nuovo Cimento 115 B (2000) 1083. 

[8] Sonnino G. Nuovo Cimento 118 B (2003) 1155. 

[9] Sonnino G. and Peeters Ph. Chaos 14 (2004) 910. 

[10] Prigogine I. 1954 Thermodynamics of Irreversible processes, (John Wiley & 
Sons) . 

[11] Glansdorff P. and Prigogine I. 1971 Thermodynamics of Structures, Stability 
and Fluctuations, (John Wiley & Sons). 

[12] Katok A. and Hasselblatt B. 1995 Introduction to the Modern Theory of Dy- 
namical Systems, (Cambridge University Press). 

[13] Fermi F.Atti R. Accad. Lincei, CI. Sc. Fis. Mat. Nat. 31 (1922) 21. 

[14] Weinberg S. 1972 Gravitation and Cosmology: Principles and Applications of 
the General Theory of Relativity, (John Wiley & Sons). 

[15] Frankcl T. 1997 The Geometry of Physics - An Introduction, (Cambridge 
University Press). 

[16] Synge J.L. and Schild A. 1949 Tensor Calculus, (Dover Publications, Inc. 
New York). 

[17] Eisenhart L.P. 1997 Riemannian Geometry, (Eight printing, for Princeton 
Landmarks in Mathematics and Physics series - Princeton University Press) . 

[18] Lovelock D. and Rund H. 1989 Tensors, Differential Forms and Variational 
Principles, (Dover Publications, Inc., New York). 

[19] York J.W. Jr. Found. Phys., 16, (1986) 249. 

[20] Horsthemke H. and Lefever R. 1984 Noise-Induced Transitions. Theory and 
Applications in Physics, Chemistry and Biology, (Springer- Verlag). 

[21] Green M.B., Schwarz J.H. and Witten E. 1987 Superstring Theory. Volume 
1 Introduction, (Cambridge University Press). 



32 



